A block Hankel generalized confluent Vandermonde 

matrix 



Vandermonde matrices are well known. They have a number of in- 
teresting properties and play a role in (Lagrange) interpolation problems, 
partial fraction expansions, and finding solutions to linear ordinary dif- 
ferential equations, to mention just a few applications. Usually, one takes 
these matrices square, q x q say, in which case the i-th column is given by 
u{zi), where we write u(z) = (l,z, ... , z' 1-1 ) 1 ^ . If all the Zi (i = 1, . . . , q) 
are different, the Vandermonde matrix is non-singular, otherwise not. 
The latter case obviously takes place when all z% are the same, z say, 
in which case one could speak of a confluent Vandermonde matrix. Non- 
singularity is obtained if one considers the matrix V(z) whose i-th column 
(i = 1, . . . , q) is given by the (i — l)-th derivative u' 1-1 ^ (z) T . 

We will consider generalizations of the confluent Vandermonde matrix 
V(z) by considering matrices obtained by using as building blocks the ma- 
trices M(z) — u(z)w(z), with u(z) as above and w(z) — (1, z, . . . , z r ~ ), 
together with its derivatives M^ k '{z). Specifically, we will look at ma- 
trices whose ij-th block is given by M' i+j '(z), where the indices i,j by 
convention have initial value zero. These in general non-square matrices 
exhibit a block-Hankel structure. We will answer a number of elementary 
questions for this matrix. What is the rank? What is the null-space? Can 
the latter be parametrized in a simple way? Does it depend on 2? What 
are left or right inverses? It turns out that answers can be obtained by 
factorizing the matrix into a product of other matrix polynomials hav- 
ing a simple structure. The answers depend on the size of the matrix 
M{z) and the number of derivatives (z) that is involved. The results 
are obtained by mostly elementary methods, no specific knowledge of the 
theory of matrix polynomials is needed. 
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Abstract 



1 Introduction and notations 



1.1 Motivation 

Vandermonde matrices are well known. They have a number of interesting 
properties and play a role in (Lagrange) interpolation problems, partial fraction 
expansions, and finding solutions to linear ordinary differential equations, to 
mention just a few applications. Usually, one takes these matrices square, qy. q 
say, in which case the i-th column is given by tt(zj), where we write u(z) = 
(1, z, . . . , z q ~ 1 ) T . If all the Zi (i = l,...,q) are different, the Vandermonde 
matrix is non-singular, otherwise not. The latter case obviously takes place 
when all Zi are the same, z say, in which case one could speak of a confluent 
Vandermonde matrix. Non-singularity is obtained if one considers the matrix 
V(z) whose i-th column (i — 1, ...,<?) is given by the (i — l)-th derivative 
u (i-i)(^)T ) or by _ iy In the i atter cas(3j one has det(V(z)) = 1. 

A slightly more general situation is obtained, when one considers the in general 
non-square q x v matrix V(z), with i-th column u^~^ (z) T , i = 1, . . . ,v — 1. 
In this case one has that V{z) has rank equal to minjg,^} and for v > q, the 
kernel of V(z) is the {y — g)-dimensional subspace of consisting of the vectors 
whose first q elements are equal to zero. Note that the building elements of the 
matrix V are the column vectors Ui(z) and a number of its derivatives. 

The observations above can be generalized in many directions. In the present 
paper we opt for one of them, in which we will consider generalizations of the 
confluent Vandermonde matrix V(z) by considering matrices obtained by using 
as building blocks the matrices M(z) = u(z)w(z) € M. qxr , with u(z) as above 
and w(z) = (1, z, . . . , z r_1 ), together with its derivatives M^ k > with < k < 
u-1. Note that M(z) = V(z) if r = 1. 

A special case of what follows is obtained by considering the matrix 

M°{z) = (M(z),...,M^ v - 1 \z)) eM. qxvr . 

In a recent paper [S], the kernel of the matrix A4°(z) (or rather, a matrix 
obtained by a permutation of the columns of M°(z)) has been studied for the 
case q = r and v < q. More precisely, the two aims of the cited paper were to 
find a right inverse and a parametrization of the kernel of the matrix M°(z). 
The two cases v — q or v < q have been analyzed in detail. For these two cases 
two algorithms have been proposed to construct the kernel and the rank and 
the dimension of the kernel have been computed. This was of relevance for the 
characterization of a matrix polynomial equation having non-unique solutions. 
The origin of this problem was of a statistical nature and lied in properties of 
the asymptotic Fisher information matrix for estimating the parameters of an 
ARMAX process and a related Stein equation. 

The results on the parametrization of the kernel obtained in [5] turned out 
to be unnecessary complicated, partly due an apparently irrelevant distinction 
between different cases, and not transparent. In the present paper we reconsider 
the problem by embedding it into a much more general approach to obtaining 
properties of matrix polynomials that can be viewed as generalizations of a 
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confluent Vandermonde matrix (in a single variable) as indicated above. The 
solution, which is shown to exhibit a very simple and elegant structure, to 
the original problem of describing the kernel under consideration, now follows 
as a byproduct of the current analysis. These are special cases of the unified 
situations of Corollary 12.51 and Proposition 14.11 of the present paper. In this 
way we generalize in the present paper the results of [8], Another approach to 
find a basis for the kernel has been followed in [7], which is also subsumed by 
Proposition 14.11 The results of [7] are closer in spirit to those of the present 
paper than the results in [5] . Both cited papers contain some examples of right 
inverses of Ai°(z), that are special cases of what will be obtained in the present 
paper. Moreover, we will generalize the results for Ai°(z) to results for the 
matrix M{z) that is defined by 

/ M{z) ■■■ M^-^(z) 
M(z) = : : 

whose blocks we denote M{z) l i , i = 0, . . . , [i— 1, j = 0, . . . , v — 1, so M.{z) 1 ^ — 
M^ + ^'{z). We will also consider the related matrix N{z) whose blocks are 

A special case occurs for the choices of the parameters r = 1, /x = 1 and 
v = q. Then we write M(z) = U q (z) and J\f{z) = U q (z), ordinary (normalized) 
q x q confluent Vandermonde matrices. The results in the next sections (almost) 
reduce to trivialities for ordinary confluent Vandermonde matrices, so the con- 
tribution of the present paper originates with allowing the parameters /i, v, q 
and r to be arbitrary. 

Classical Vandermonde matrices and confluent Vandermonde matrices have 
often been studied in the literature, see [9j|3] for definitions. More recent papers 
often focussed on finding formulas for their inverses and on efficient numerical 
procedures to compute them, a somewhat random choice of references are the 
papers [TJ |2j HI HI E [10l HU HI] ■ Nevertheless it seems that the generalization 
A4(z) of the confluent Vandermonde matrix is, to the best of our knowledge, 
unknown in the literature. 



1.2 Notations and conventions 

Derivatives of a function z H> f(z) (often matrix valued) are denoted by f( k \z) 
or by (g^) f{ z )- The variable z is in principle C-valued. If A is an m x n 
matrix, for notational convenience we adopt the convention to label its elements 
Aij with i — 0, . . . , m — 1 and j = 0, . . . , n — 1. We will see shortly why 
this a convenient convention. Entries of a matrix are also indicated both by 
superindices, according to what is typographically most appropriate. When 
dealing with matrix polynomials A{z), in proving results we usually suppress 
the dependence on z and simply write A. 

If A is a block matrix, we sometimes use subindexes to indicate its consti- 
tuting blocks, but more often denote its blocks by superindices, so we write A lJ 
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or A 13 in the latter case. Also single supcrindices are used and we will come 
across notations like A k . These should not be confused with powers of A. The 
meaning of A k will always be clear from the context. 

Throughout the paper q, r, /i and v arc fixed positive integers, although often 
certain relations among them are supposed (e.g. q + l<v<q + r). Given 
the integers q, r, fj, and v and the variable z, we consider a number of matrix 
polynomials. The first are u(z) = (1, . . . , z 9_1 ) T and w(z) = (1, . . . , z r_1 ). 
With the above labelling convention, valid throughout the paper, we have for 
the elements of u{z) the expression Ui(z) — z % for i = 0, . . . , q — 1 . Next to these 
we consider the matrix polynomials 

M(z)=u(z)w(z) el» xr , 

M%z) = M^(z) = (-^-yM(z) e R qxr , 

Ai (z) = (^(z),...,^_ 1 (z))eR^, 

Jtf{z) = \m^{z) = -(^-yM{z) G R qxr , 
]■ .]■ dz 

N\z) = «(*), . . . , A/?-! (z)) e R^. 

We also consider the [iq x vr matrix polynomial M(z) which is defined by its 
zj-blocks 

M ij (z) =MW>(z), (1.1) 

for i = 0, . . . n — 1, j = 0, . . . , v — 1. The size of A4(z) is equal to \xq x vr. If 
fi = 1, we retrieve M°(z). Note that M(z) has a block Hankel structure, even 
for the case [i^v. 

Along with M(z) we consider the matrix M{z) that has a block structure 
with blocks 

for i = 0, . . . , [i — 1 and j = 0, . . . v — 1. Note that AA(z) also has dimensions 
/zg x vr and that A/"(z) reduces to J\f°(z) for /i = 1. Unlike Al(z), A/"(z) doesn't 
exhibit a block Hankel structure. The following obvious relation holds. 

M(z) = (Dp <8) I q )N(z){D v ® I r ), (1.2) 

where D M is the diagonal matrix with entries Dj, = i!, i = 0,...,/z — 1 and D v 
likewise. 

Let us now introduce the matrices U q (z) (of size q x q) and W r {z) (of size 
r x r) by 

u,(«) = («(«),..., u^- 1 )^)) 

and 

VK r (z) T = (w(z) T ,...,w (r - 1) (z) T ). 
Along with the matrices U q {z) and W r (z), we introduce the matrices 

U q (z) = (u (z),...,u q - 1 (z)), (1.3) 
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with 



U j( Z ) = — U 3 = °) •••)?- 



and Wr(,z) given by 



Wr(«) 



T 



(wo(z), • ■ • , w r _i(z)), 



(1.4) 



with 



= — w {3 '(z), j = 0, . . . ,r. 



We have the obvious relations 



U q {z) = U q (z)D q , W r (z) = D r W r (z), 



where D q and D r are diagonal matrices, similarly defined as D^. One easily 



that Uglz)" 1 = U q {—z). Similar properties hold for W r (z). Later on we will 
also come across the matrices U^z) and W v {z), which have the same structure 
as U q (z) and W r (z) and only differ in size (unless /i = q and v = r). 

Finally we introduce the often used square shift matrices Sk £ R fcxfc (fc 
arbitrary), defined by its ij-elements Si+ij (Kronecker deltas). Other matrices 
will be introduced along the way. 

The structure of the remainder of the paper is as follows. After having fixed 
some notation and other conventions, in Section[2]we shall introduce an auxiliary 
matrix polynomial A{z) that is instrumental in deriving properties of A4(z) 
and J\f(z). Among them are factorizations, which will be treated in Section |3J 
These lead to establishing the rank of A4(z), which is shown to be independent 
of z. A major issue in the present paper is to find a simple and transparent 
parametrization of the kernels of the matrices A4(z) and N(z). This is done 
first in Section |4] for the special case of J\f° (z). The results of that section will 
be used in Section [5] to characterize the kernel of Af(z). The results of Section [6] 
form the cornerstone of finding right inverses of M(z), especially those that are 
not dependent on z, which is a main topic of the final Section [7] 



This section is devoted to the matrix A(z), to be defined below, that is in- 
strumental in deriving properties of A4(z) and N{z). In particular it is used 
for obtaining useful factorizations in Section [3] The matrix A(z), also of block 
Hankel type and of size nq x vr is defined by its blocks Aij(z) :— A l+3 (z) (here 
i + j is used as a super index) of size gxr, for i = 0, . . . /J, — 1, j = 0, . . . , v — 1. 
The matrices A k (z) for k > are the q x r quasi-symmetric matrix polynomi- 
als (this matrix polynomial is only square for q = r) given by their elements 



verifies that U q {z) has elements U q 3 (z) — z % 3 



(*.) and by a simple computation 



2 The matrix A(z) 



= 0,...,q- 1) 




(2.1) 
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Note that fc in A k {z) is used as a super index and in z k as a power. The block 
Hankel matrix A(z) of size /j,q x vr is then defined by its blocks A^(z) which 
are equal to A t+: > (z), for i = 0, . . . , fx — 1 and j = 0, . . . , v — 1. One easily verifies 
that (with A^-i. = = 0) for k > 1 it holds that 

4(z) = zA\-\z) + A^liz) + A k -\{z). 

In matrix notation this relation becomes 

A k (z) = zA k - 1 {z) + A k ~ 1 {z)S r + Sj 'A k -\z). (2.2) 

Let J q (z) = zl q + Sj . Then we can rewrite (|2.2[) as 

A k (z) = J q {z)A k - x {z) + A k - 1 {z)S r , (2.3) 

or as 

A k {z) = SjA k - 1 (z) + A k - 1 (z)J r (z) T . (2.4) 

By induction one easily proves that the recursion (|2.3I) leads to the explicit 
expressions, involving powers of Sj and S r , 

A k {z) = j2( k )w i A°s k r 



and 

fc 



A k +\z) = £(*)j,(*)*A l (*)S*- 



(2.5) 



Likewise one shows that (|2.4I) leads to 

A k (z) = ( k ) (SjyA (J r (zy) k -\ (2.6) 



On the other hand, the recursion (|2.2I) has solution 



Ak ^ =p o ( k j) zk ~ 3 i (^(sjr^sr- (2.7) 

It also follows that 

A k (0) = ( k ) {SjYA°S k -\ 

i=Q ^ l ' 

Look back at the A k {z) as defined in (I2.1j) . One computes for fc > 

Henceforth we shall write A k instead of A k (0). Note that A k , which is in general 
not square, has only nonzero elements on one 'anti-diagonal', the one with ij- 
elements having sum i + j = k. We have the following result. 
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Proposition 2.1 The matrices A k = j4 fe (0) satisfy the recursion for k > 



A k+i = A k S r + SjA k , 



and hence 

k 



A k = J2(^J(SjyA°S*- j . 



3=0 

Morever, we also have the polynomial expansion 

k 

.\z k - j A j . 

3=0 ' ' 



Proof The recursion follows by taking z = in (|2.2I) . The expansion follows 
from Proposition ^. II and Equation (|2.7p . □ 



Two more matrices will be introduced next. First, let A be the [iq x vr block- 
matrix whose ij-th block (of size q x r) is given by A 1 ^ = (Sj) J A S* for i — 
0, . . . , \i — 1 and j = 0, . . . , v — 1. 

Second, we introduce the matrix C^^ (z) of size y,q x /iq with blocks C fl . q (z)ij = 
{j)Jq( z Y ~^ f° r * — J ana - zer0 e ^ se ' = 0,...,/i — 1. Since is block 

lower diagonal with identity matrix as diagonal blocks, it follows that C flljq (z) is 
invertible and that its inverse has ij-block equal to (1) Jq(zy~ J (— Likewise 
one defines the matrix C u . r {z) of size vr x vr. 

Theorem 2.2 TTie factorization 

A(z) = C^ q {z)AC y , r {z) T (2.9) 

raoZds true. For z — and ^4 = .4(0) and £ M , g = £ Mi9 (0) and £„, r = £^,.(0), 
i/iis becomes 

A = C^ q ACl r . (2.10) 

Moreover, det(,4(z)) = det(.A), when both matrices are square. 

Proof We compute the r/-block on the right hand side of (|2.9[) . Using the 
definitions of C fJljq (z), C Ujr (z) and A, we obtain (see the explanation below) 

(£ M , 9 (z)^C, r (z) T ) y . - J2 C ^ z ^ kl £„Az)7j 

k,l 



iilQ^y-'isjy^s^^iJAzfy- 1 

k=0 1=0 ^ ' ^ ' 

E Q)^w i -*(E(5, T ) , A°(j)(j P (z) T y- l )5 r fc 

fe=0 ^ ' Z=0 ^ ' 

j2(l) j ^y- k A i (z)s k 

l — n V / 



k=0 

= A l+: >(z). 
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In the third equality above we used that J r (z) T and S r commute, in the fourth 
equality we used (|2.6p with an appropriate change of the indices, whereas the 
last equality similarly follows from (|2.5p . The relation between the determinants 
follows from det(£ Mig (z)) = 1. □ 

Let e, be the i-th standard column basis vector of M. q , i = 0, . . . , q — 1 and let /, 
be the i-th standard column basis vector of W, i = 0, . . . , r — 1. For convenience 
of notation, we put = for i > q and fj — for j > r. With this convention, 
we always have for example (SjYea = e^. 

Theorem 2.3 It holds that the block A tJ — CjfJ , the rank of A is equal to 
min{/x,r} x min{^, q}. If A is a fat matrix, i.e. \iq < vr , then A has full 
(row)rank iff fx < r and v > q. If A is a tall matrix, so y,q > vr, then it has 
full (column) rank iff fj, < r and v < q. In the special case that A is square, so 
fiq = vr, we get that A has full rank, and it is then invertible, iff fj, = r and 
v = q, in which case A is a matrix of size qr x qr. In this case we have for the 
inverse (A)^ 1 — (A) T and det(„4) = (— l)4<? r (9- 1 )( r - 1 ). Specializing even more 
toq = r, we get det(A) = (-l)^- 1 ) . 

Proof Since A" = e / T and A ij = {Sj) j A°S'i., it follows that A lj = fjej. 

The rank of A is equal to the rank of A(A) T , which is easy to compute. We 
get for its ij-th block J2i=o Zlf7 fj e J = fi fj ELo e i e 7 ='■ fi 'fj^v-i, where 
Iq.v-i = ^}2d = o e i e l ■ ^ follows that rank(/ gj „_i) = minjg, v}. Furthermore we 
have fj fj — for i ^ j and fj fi = 1 iff i < r — 1. We conclude that A(A) T is 
block diagonal, where the diagonal M-blocks are equal to I q , v -i for i < r — 1 and 
zero otherwise. The number of nonzero diagonal blocks is equal to min{/x,r}, 
hence the rank of -4(,4) T is equal to min{/i, r} x min{z/, q}. 

Assume that A is fat and that \x < r and v > q. Then the rank of A 
equals fiq, the number of rows of A. For the converse statement we assume that 
fi > r (the case v < q can be treated similarly). In this case the rank becomes 
r x minji/, q}, which is strictly less than fiq, the number of rows of A, which 
then has rank deficiency. The dual statements for a tall A follow by symmetry. 

Assume next that the matrix A is square. It has full rank iff the two sets of 
conditions for the tall and fat case hold simultaneously, which yields the asser- 
tion on invertibility. Assume then that fi — r and v = q. By the computations 
in the first part of the proof we see that the diagonal blocks of _4(.4) T are all 
equal to I q . Since there are now r of these blocks, we obtain that ^4(^4) T is the 
qr x qr identity matrix. 

To compute the determinant for this case, we observe that the columns of A 
consist of all the basis vectors of M. qr , but in a permuted order. Therefore, the 
determinant is equal to plus or minus one. To establish the value of the sign, we 
compute the number of inversions of the permutation. This turns out to be equal 
to \qr{q — l)(r — 1), which results in a determinant equal to +1 iff this number 
is even, and —1 in the other case. If r = q, then (— l)!9 2 (<?-i)" _ j± 
way of computing the number of inversions is to write the order of the numbering 
of the column basis vectors in rectangular array. Decomposing every number 
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x 6 {0, 1, . . . , qr — 1} in a unique way as x = nq + m with m £ {0, . . . , q — 1} 
and n 6 {0, . . . , r — 1}, we can identify every such x with a pair (m, n). An 
inversion i(x,y) occurs when x < y, but in the permuted order x is preceded 
by y. For every x the number of inversions i(x,y) is the number of elements 
in the rectangle strictly to the South- West of x in the rectangular array. So, 
if x corresponds to (m, n) then the number of inversions i(x, y) is equal to 
(q — 1 — m)n. Summing these numbers for m = 0, . . . , q — 1 and n = 0, . . . , r — 1 
yields the total number of inversions. □ 

Remark 2.4 In the case were A is square and invertible, it is the permutation 
matrix having the property vec(X T ) = Avec(X), for any X S R 9Xr . 

Corollary 2.5 The rank of the matrix A{z) is for all z € C egttaZ £o min{/Lt, r} x 
min{z/, 5}. If [iq = vr, then A{z) is square and invertible and det(A(z)) — 

(_l)l9»-(9-l)(r-l) > 

Proof The assertion on the rank follows from Thcorems l2.2l and l2.3l upon noting 
that the matrices C f±tq (z) and C v>r {z) are invertible. Since det(C l±t g(z)) — 1, also 
the assertion about the determinant follows. □ 

3 Factorizations of the matrices A4(z) and N(z) 

In the present section we obtain a factorization of the matrix polynomial A4(z), 
from which a factorization of the matrix Af(z) follows as a simple corollary. In 
the next proposition we use the Kronecker symbol (g> to denote tensor products. 

Proposition 3.1 The factorization 



holds true. Moreover, M.[z) and A have the same rank equal to min{/i, r} x 
min{t/, r}. 

Proof Computation of the product (I q (g> U q (z))A(I q <8> W r (z)) by using the 
definition of A as block matrix, yields a matrix that consists of blocks 



M(z) = (J M ® 17,(«))w4.(J„ <S> W r («)) 



(3.1) 



g— 1 1 — 1 



U q (z)A k W r (z) = 



i=0 j=0 



Using Equation (|2.8p . we get that this expression reduces to 



9-1 r ~ 1 (g-l)Afc - . 

EE r^t^w^-i = E (Jjuww^w 
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since u^(z) = for i > q and w^ k ^(z) — for i < k — r. Recall that 
the blocks M.(z) consist of the derivatives 4- l; {u{z)w{z)). The product rule for 
differentiation then yields (|3.ip . The statement concerning the rank immediately 
follows from Theorem 12.31 since U q (z) and W r (z) are invertible. □ 

Remark 3.2 The factorization in Proposition 13. 1 1 exhibits a nice form of sym- 
metry. Therefore it would be nice if also the matrix A could be factorized in 
some symmetric way. There doesn't seem to be an easy way to do this. Con- 
sider a truly symmetric case, for instance [i = v = q = r = 1. In this case A is 
non-singular and we have 



A 



^10 1^ 
10 
10 

\1 2/ 



A first reasonable factorization in this symmetric case would be of the form 
A = AA T , which would imply that A is positive definite. One easily sees that 
this is not the case. As a next attempt, one could try to use the singular value 
decomposition of A, but the eigenvalues of A are not particularly nice, so this 
look as a dead end too. 

Remark 3.3 Instead of the matrices U q (z) and W r (z), one can also use the 
matrices U q (z) (see (|1.3p ) and W r (z) (see (|1.4[0 to get a factorization of M.(z). 
One then has to replace the matrix A with A, whose blocks Aij are equal to 
where the matrices A k are specified by their elements 

~ fe ( fc! if % + j = k 
v I else. 



M{z) = (/„ <g> U q {z))A{I v ® W r (z)), (3.2) 



One then gets 



which can be proved in the same way as (|3.1[) . or by using this identity and the 
relation 

A k = D q A k D r . (3.3) 

Next we derive a factorization of the matrix J\f(z). Let the matrix A be defined 
by its /i x v blocks A^ , for i = 0, . . . /i — 1, j = 0, . . . , v — 1, where the matrices 
$3 e Rixr have elements 

A i j= ( C+ J ) if k + l = i + j 
kl \ else, 

for k = 0,...,q- 1, I = 0,...,r- 1. 
Proposition 3.4 The factorization 

N{z) = {I„ ® U q {z))A{I v ® W r {z)) (3.4) 
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holds true. The matrices A and A are related through 

(D„ ® I q )A{D v ® I r ) = (I„ ® Dq)A(I v ® D P ). (3.5) 

Further relations between A, A and A are 

A = {Dfj, <8 / g )^(-D„ (8 Ir) (3-6) 
= (!„ ® A,).A(I V ® A-) (3.7) 

Proof The proof of the factorization p. 41) is similar to the proof of (|3.1[) . Re- 
lation (I3.5[) follows by an elementary computation. Likewise one proves (|3. 61) 
and JO- □ 



Proposition 3.5 Additionally we have the following factorizations of A4(z) 
and M(z). 



A{z) = 


{U^{z) ® I q )A{Q){W v {z) ® I r ) 




(3.8) 


M{z) = 


(I^®U q (z))M(0)(I„®W r (z)) 




(3.9) 


M{z) = 


(%{-z) ® U q (z))A(z)(W„(-z) 8 




(3.10) 


N{z) = 


(U^z)- 1 ®U q {zj)A{z)(W v {z)- x 


®W r {z)) 


(3.11) 


Af(z) = 


(I f t®U q (z))N(0)(I u ®W r (z)) 




(3.12) 



Proof First we show an auxiliary result. Let V^ q {z) = C IMtq (z)C IJ , iq (0) 1 . Then 

V M (z) = U li (z)®I q . (3.13) 

Indeed, by direct computation starting from the definition of C^ tq (z), one finds 
that the ij-block V^. q (z) %3 equals ^z % ~^l q for < j < i < fi — 1. Since 

Ufj,(z)ij — (y)z % ~ 3 , the result follows. 

We now proceed to prove the indentities. Use (|2.9p . both for arbitrary z and 
z = 0, and Equation (|3.13p to get (|3.8j) . To obtain (|3.9j) . in the same vain we 
use (|3.1|) twice. Then (|3.10j) follows upon combining (|3.8|) and (|3.1|l . Finally, 
(f3~TT|) follows from (f3~T0f and JO) and JSH]) follows from JO and JO]) . □ 

In the situation where .M(z) and N(z) are square, we are able to compute their 
determinants. 

Corollary 3.6 Let fiq = vr. Then all relevant matrices are square and we have 
the following expressions for the determinants of A4(z) and Af(z). 

q—l i — 1 

det(M(z) = (JJi!)" (Y[j\)" det(A), 

det(7V(z)) = det(i), 

ii— l i/— l 
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with det(„4) given in Corollary \2.5\ Hence the matrices M(z) and J\f(z) are 
unimodular. 

Proof The matrix Iy,®U q (z) has determinant one, which gives the first equation 
in view of p.4[) . Use (I3.2[) . (|3.3[) and its consequence 



det(.A) = (det(.D g ))^(det(I> P )) , 'det(.A) 
to get the formula for det(Ad(z)). Then the last equation follows from (|3.5p . □ 



In the next sections we investigate and characterize the kernel of the matrix 
Af(z) for different values of the parameters. It turns out that characterizing 
the kernel of Af(z) yields more elegant results than characterizing the kernel 
of M(z). On the other hand, a factorization of A4(z) leads to more elegant 
expressions than factorizations of Af(z). Of course, in view of (|1.2j) . results for 
one of the two can easily be transformed into results for the other. We will 
focus on the kernel of Af(z) only, simple because the obtained results have a 
more elegant appearance. As it turns out, it is instrumental to consider the 
kernel of Af°(z) (which is Af(z) for /x = 1) first, since the results obtained for 
this case, serve as building blocks for the kernel of M{z). 

4 Characterization of the kernel of J\f°(z) 

We will investigate the kernel of J\f°(z) for different values of v, q and r. This 
is in different notation the problem alluded to in the introduction and earlier 
investigated in [5] for a special case. Notice that Af°(z) is q x rz/-dimensional. 
Furthermore for v < q all derivatives u^ k '(z) (k < v — 1) are nonzero and 
moreover linearly independent vectors, whereas vS k \z) = for k > q. It follows 
that the set span{u(z), . . . , u u ~ 1 (z)} has dimension equal to min{z/, q}. It is then 
an easy exercise to directly see that the rank of Af°(z) is also equal to min{j/, q}. 
Of course, this also follows from Proposition 13.41 In order to characterize the 
kernel of Af°(z), we have to discern three different cases. These are v < q. 
q+l < v < q + r and v > q + r, each case being treated in a separate section. In 
all cases, our aim is to find simple, transparant paramctrizations of the kernels. 

4.1 The case v < q 

Let F(z) be the r x (r — 1) matrix given by 



l-z 
1 -z 















1 



F(z) 







— z 



\0 



1/ 
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Observe the trivial but crucial property, that the columns of F(z) span the 
(r — l)-dimcnsional null space of w(z). 

Next we consider the matrix K(z) that has the following block structure. It 
consists of matrices (z) (i = 0, . . . , v — 1, j = 0, . . . , v — 1) where each (z) 
has size r x (r — 1) and is given by 



K u {z) 
K iii+ i(z) 



F(z), 
-F'(z), 



whereas all the other blocks are equal to zero. For r = 1 the matrix K(z) is 
empty and by convention we say that the columns of K (z) in this case span the 
vector space {0}. Note that K{z) has dimensions vr x v(r — 1). 

The matrix K(z) looks as follows, where we suppress the dependence on z and 
omit zero blocks. 



K = I V ®F + S„®F' = 



F' 
F 



F' 



V 



F' 
F) 



(4.1) 



where I v is the ^-dimensional unit matrix and S v the ^-dimensional shift matrix, 
with elements = S i+ ij. 

Proposition 4.1 Let v < q. The vr x v(r — 1) matrix K{z) has rank v(r — 1) 
and is such that AT°(z)K(z) = 0. In other words, the columns of K(z) form a 
basis for the kernel ofAf (z). 

Proof Pick the j-th block column of K (z), (z) say (j e {0,...,v— 1}). We 
compute N Q {z)Ki{z). For j = 0, this reduces to M(z)F(z), which is zero, since 
w(z)F(z) = 0. Then, for j > 1 we have Kj_ij(z) — F'(z). Hence, for j > 1 we 
get (using that F^iz) = for k > 1, 



Af°(z)Ki(z) = J2Mi(z) K ii(z) 



i=0 



= -M^{z)F{z) + - 1 - M^- 1 \z)F'{ 



(.1 1)! 



-(M {l \z)F{z) + jM^- 1 \z)F'(z)) 



-(M(z)F(z)) 



U) 



= 0. 



Hence Af°(z)K(z) — 0, so all columns of K(z) belong to kcr7V°(2;). Since 
we know that the rank of A/" (z) is equal to v, we get that dimker Af°(z) = 
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rv — v = (r — l)v , which equals the number of columns of K(z). Since K(z) is 
upper triangular with the full rank matrices F(z) on the block-diagonal, it has 
full rank. Therefore, the columns of K(z) exactly span the null space of TV (z). 

□ 

Remark 4.2 A nice feature of the matrix K(z) in Proposition 14.11 is that it 
is a matrix polynomial of degree 1, only whereas the matrix polynomial J\f (z) 
has degree q + r — 2. 



4.2 The case q + l<v<q + r 

Next we consider what happens if q + 1 < v < q + r — 1. Note that this case is 
void if r = 1. Henceforth we assume that r > 2. Let the matrix G(z) £ M. rxr 
have elements Gij(z) = z^ % ~ x for j > i and zero else. With S — S r the r- 
dimensional shift matrix, we have the compact expression G(z) = S ^2^L (z S) k , 
since S k = for k > r. We now give some auxiliary results. 

Lemma 4.3 Let G(z) be as above and let S be the shifted r xr identity matrix, 
Sij = Si + ij. It holds that 

G{z) = S(I - zS)- 1 (4.2) 
-,G^{z) =S j+1 (I -zS)-!' 1 . (4.3) 

r- 

-G^ [z) = G(zy +1 (4.4) 



— L- G^ +1 Hz) = -G^{z)G(z) 1 (4.5) 

u + i)! r- 

and the matrices G^ (z) and G?W (z) commute for all i,j > 0. One also has the 
following two equivalent properties 

-w {3 \z)=w{z)G{zy . (4.6) 

r- 

-w {3+l \z) = -w {3) {z)G{z) 1 . (4.7) 



(J W j\ 



Moreover, 



(w(z)G(z))^ = (k + 1)! w(z)G(z) k+1 . (4.8) 
Proof To prove (|4.2I) . we recall that = for j > r. Hence 

oo 

G(z) = S^2(zS) k = S(I - zS)-\ 

Then (|4.3j) simply follows by differentiation, which immediately yields (|4.4[) , 
from which (14.51) also follows. 
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We next prove (|4.6[) for j = 1. This is then equivalent to w'(z)(I — zS) = 
w(z)S, which can be verified by elementary calculations. We proceed by induc- 
tion and assume that (I4.6[) holds. Then we have 

' w^ +1 \z) = -^—-^w^(z) 



+ v ' (j + l)!dz 

j! d 



(w(z)G(zY) 



(j + 1)! dz 

-(w'(z)G(zy + ] w{z)G{zy- 1 G'{z)) 



(.1 + 1)! 

+ jw(z)G(z) j+1 ) 

3+1 



= w(z)G(z) 

This proves (|4.6p . which is easily seen to be equivalent to (|4.7[) . Equation (|4.8p 
follows by induction, using (I4.6|) . □ 

Remark 4.4 In fact, above statements about G(z) also follow from the parallel 
ones concerning w(z) and the relation 

(w(z)S 
; 
w(.z)S r 

Indeed, since the matrices G(z) and S commute, one has for instance 

1 

j! 

which together with similar relations leads to (|4.4I) . 
Lemma 4.5 One /las for j > 1 

i(M«(«) - u (j) (z)w(z)) = _L_A/(J- 1 )(z)G(z). (4.9) 

(j - 1)! 



-«/%)S* = w(z)SG(z) i , 

7 I 



For j > i/iis reduces to 

jMU\z) = —^M^\z)G{z). (4.10) 

More generally, one has for m,p > 

M (z) = M ^Z G ( Z )P = ?! ^2 G(z )f>+™ (4.11) 

(q-l+p + m)\ (q-l + m)\ w 
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Proof Recall that M(z) = u(z)w(z). Below we use (|4.7[) to compute 



^M {j \z) = \Y1 ( ] A )u {:J ~ l) {z)w {t) {z) 



v. * — ' V l 

J j=0 



2^ a 



(j-i) 



G(z) 



-ui(z) 



(j 


-1)! 




1 


U 


-1)! 




1 


U 


-1)! 



?E( ] l _ l M J - l \z)w^- 1 \z)G(z) + 



w(z) 



3-1 



1=0 



1 M^ 1 \z)G{z) + ^-P- W (z). 



Equation (|4.10p follows from (|4.9j) . because = for j > q. Finally, (|4.11[) 
follows by iteration of (|4. 10[) . □ 

We consider the matrix polynomial K{z) which for the present case has the 
following structure. For i, j = 0, . . . , q — 1 it has blocks Kij(z) of size r X (r — 1) 
such that Kjj(z) — F{z) and Kj-x,j(z) = F'{ z ), F(z) as before. For i, j = 
q, . . . , v — 1 the blocks are Kjj(z) = I r and Kj-\ j(z) = —G(z), all of size r x r. 
Finally, we have that K q -\ iq {z) = —G(z). All other blocks are equal to zero. 
Notice that K(z) is of size vr x (yr — q). One easily verifies that K(z) has full 
column rank. 

The matrix K(z) looks as follows, where again we suppress the dependence on 
z and omit zero blocks. 



( F F' 

F F' 



K 



F' 
F 



\ 



-G 

/, G 

L -G 



-G 

Ir ) 



(4.12) 
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A compact expression for K is 

R = fl q ® F + S q ® F' -i q ij- q ® G 

V ® J r - <g> G 

where f„_ q is the first standard basis vector of R 1 ^ 9 and l q the last standard 
basis vector of R 9 . 

Proposition 4.6 Le£ g + l<i'<g + r — 1. The vr x (Vr — q) matrix K(z) 
is such that Af° (z) K (z) = 0. In other words, the columns of K(z) form a basis 
for the kernel ofJ\f°(z). 

Proof Next we proceed as in the proof of Proposition ^. II We pick the j-block 
column K j {z). If 0_< j < q - 1, then Kjj(z) = F{z) and Kj-ij(z) = F'(z), 
whereas all other Kij(z) are zero. The computation of Af°(z)K 3 '(z) is then 
exactly as in the previous proof. Next we consider the case where q < j < v — 1, 
which is quite different. Again we pick the j-th block column of K(z). Recall 
the definition of the K^{z) for this case. We get, 

N\z)Ki{z)=Y J ^(z)K l3 {z) 

= ^M^(z)-j-±-^M^- 1 Hz)G(z) 
= 0, 

in view of Equation (|4.10l) . This shows that K(z) belongs to the kernel of Af°(z). 
Since J\f°(z) has rank q, the dimension of the kernel is equal to vr — q, which is 
equal to the rank of K(z). Hence the columns of K(z) span this kernel. □ 

Post-multiplying the matrix K(z) by 



/ ir-l 






\ 


Ir-l 


■• 

Ir-l 










I r G ■■■ ■ 


■ G u - q - 1 






I r G ■ 


■ G v - q - 2 


V 






G 

Ir J 



(4.13) 
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we obtain 



{ F F 1 






\ 


F F' 


■. F' 








F 


-G -G 2 ■■■ 


■■■ -G v ~ q 






Ir 

Ir 




\ 






Ir 1 



an alternative matrix whose columns span kerA/"°(z). 
4.3 The case v > q + r 

For this case, the kernel ofAf°(z) is closely related to what we have obtained in 
the previous case. 

Proposition 4.7 Let v > q + r. Consider the matrix K*(z) of {^.12^ in the 

special case that v = q + r — 1, then K*(z) has (q + r)(r — 1) columns. We have 
that kerj\f a (z), which is now vr — q- dimensional, is the product of the space 
spanned by the columns of K*(z) andW( u+1 ~ q ~ r \ 

Proof Since the highest power of z that appears in Af°(z) is q + r — 2, we have 
that the matrices M"'(z) arc identically zero if v > q + r — 1, whereas the 
matrix (^M(z), . . . , ^ q+ ^_ 2 y_ M ( - q+r ~ 2 ' ) (z)) is the same as the matrix J\f°(z) for 
the case v = q + r — 1. The assertion follows by application of Proposition 14.61 

□ 

5 Characterization of the kernel of M{z) 

We have seen that we had to distinguish three different cases to describe the 
kernel of Af°(z). The same distinction has to be made in the present section. 

5.1 The case v < q 

First we introduce some more notation. For the matrix K(z) of (14. ip we now 
write K (z) and for k = 1, . . . , r — 2, we define Kf.(z) in the same way as K(z), 
but now with F replaced with F^{z), a (r — k) x (r — 1 — k) matrix having the 
same structure as the original F(z) of Section [4. 11 so Fk(z)ij = —z, if i = j and 
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Fk{z)ij = 1, if i = j + 1. In particular Fq(z) = F(z). Formally, for 1 < j < r — 2, 
we have 

F h (z) = (l r - k {r - k)xk )F(z)( n Ir - k - 1 ), 

whereas Fk(z) is the empty matrix for k > r — 1. By K^~ l {z) for /i < r, we 
denote the product K (z)K 1 (z) ■ ■ ■ K^-i(z) of size vr x v{r — fi), whereas we 
take lC^~ l {z) the zero matrix for fx>r. 
For j = 0, . . . , r — 1 we put 

I r —j 

and Mj(z) = u(z)wj(z). With this convention, we have wq = u>, Mq = M. 



Wj (z) = w(z)[ Q ^ )=(1,z,...,z<-'-i) 



Lemma 5.1 Let w\(z) = (l,z,...,z r 2 ) and M\(z) = u(z)w\(z). For j > 1 
one has 

w (3 \z)F{z) = -jw < - J - 1) {z)F , (z) = jw¥~ 1) (z) (5.1) 

M^{z)F{z) = -jM^- 1 \z)F'{z) = jM[ j - 1 \z). (5.2) 

Proof Since w(z)F(z) = 0, also ■£ I (w(z)F(z)) = for all j > 1. It then follows 
that (z)F(z) +jw < j- 1 ){ z )F'{z) = 0, since all higher order derivatives of F 
vanish. Using — w{z)F'(z) = wi(z), we arrive at (|5.ip . Similarly, one obtains 
(1531). □ 



Theorem 5.2 The kernel of N(z) is spanned by the columns of IC^ 1 (z) and its 
dimension is equal to v(r — /i) + . Hence the rank ofAf(z) is equal to vimii{fj,,r}. 
In particular the matrix M{z) has full rank iff /x > r. 

Proof The arguments used in the proof of Proposition 14. 1 1 can also be applied 
to this more general case. Let for k > 

^ (z) = ^^ M(fc) '---'(^ M(fc+ " 1)(z) )- 

Proposition O yields M°(z)K (z) = 0, and therefore N°{z)K,^- 1 {z) = 0. 
Consider now k > 1. As before, K^(z) denotes the j-th block-column of 
K (z) = K(z). Then, from (|5~2|) it follows that 

N k {z)K*{z) = l M ^(z)F(z) = ^-L-Mf" 1 ^). 

For j > 1 we get 

N k {z)Ki{z) = l(_L_Af( fe +.?- 1 )(z)F'(z) + -KI ( - k +^(z)F(z)) 

kl {j-iy- r- 



(k-iy.j\ 1 



M{ ' {z), 
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where we used (|5.2p again. It follows that for k > 1 

Af k {z)K Q {z)=Mt\z), (5.3) 

where 

which is a matrix of size q x v(r — 1). Invoking Proposition 14. 1 1 again, we obtain 
for k = 1 

AP-WKoWK^z) = Af?(z)Ki(z) = 0, 

and hence M 1 {z)K tl ,-\{z) = 0. Assume that \i < r. To have Af(z)fC fl ~ 1 (z) = 0, 
we need Af (z)K,^ (z) = for k = 0, . . . ,fi — 1. This now follows by iteration 
of (|5.3[) . In fact, by induction, one can show 

Ar k (z)K (z) ■ --Kj-xiz) = Af*- j (z) for j < k, (5.4) 

where 

W = (^-yy (^Mf - ) , . . . , ^Mj*-'- 1 ' (*)) . 

For j > k + 1 one has M k (z)K {z) ■ ■ ■ Kj-i(z) = 0. 

Since each of the matrices Kk(z) for k < r — 1 has full rank, which is equal 
to v[r — k — 1), we get that K,^~ 1 (z) has rank equal to v[r — All assertions 
for /j, < r now follow. On the other hand, for fi > r, the matrix M{z) has rank 
equal to tr, and therefore has a zero kernel. □ 

Remark 5.3 The matrix K ti ~ 1 (z) for /i < r — 1, which is of size x v(r — \i) 
turns out to be upper block-triangular. Consider for this case the product 
— Fo(z) ■ ■ ■ F^i(z). Then one has for j > i the ij-block 

which can easily be shown by induction. 
5.2 The case q + l<u<q + r 

Next we extend the result of Proposition 14.61 to obtain the kernel of the matrix 
Af(z) for the present case. The approach that we follow is the same as the one 
leading to Theorem 15.21 

To obtain our results, we need to introduce additional notation. Let Gkj(z) 
denote the upper left block of G{z) having size (r—k) x (r—j) for < k, j < r—1. 
So 

Gkj(z) = (I r -k O(r-k)xk) G(z) I r ~ 3 ]. 

yJ 3 x(r~j)J 

We also need the matrices Qo(z) — I r , and Gi{z) = Gi^-i(z) ■ ■ ■ G\o{z) £ 
K (r-<)xr for < i < r- 1. 
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Lemma 5.4 It holds that Gjj(z)Gjj-i(z) = Gjj-i(z)Gj-ij-i(z) for j > 
1 and for j > i > k > one /ias Gji(z)Gik(z) — Gj t i-\(z)Gi-i t k(z) — 
G 3 k(z)Gkk(z)- 

Proof The first assertion follows from the decomposition 

C - ( G 3>J-A - ( G 33 9 

Gj - lj - x ~ ^ o ) - \ 

where g is the last column of Gj t j-%. For the proof of the second assertion 
we need the following property of the shift matrix S € M. k [k according to the 
context): IS = S, where 

0, 



I = 



Since any Gu{z) is of the form S(I — zS) 1 (with S of size (r — i) x (r — i), see 
Lemma l4.3p . we have IqGu(z) — Gu(z). Then we compute 



G Jt G lk = (I 0) Q (/ 0) ( 

= (7 0) Gj-i^-i ^ Gi-\ t i- 



— [I 0) Gj-i^-i/Gj-i^-i 

= Gj,i—lGi—l,k- 



□ 



Lemma 5.5 For i > I, it holds that 



-G$(z) = G i0 {z)G{z)\ 

Gi(z) = Gio^Gizy- 1 (5.5) 
G«(2)&(*) = £(z)G(z). (5.6) 

i a) 

Proof Using the definition of Gjo, the equality ^G\ q = G^G 1 immediately 
follows from (|4.4|) . The second equality (|5.5p is obviously true for i = 0. We use 
induction. Let i > 1 and assume that ft = GioG z_1 . Then, using Lemma |5.4[ 
g i+1 G = G l+ i, i g t G = G i+ i ti GioG i = G i+1 ,oG i+1 . To prove we use J52) 

and Lemma El to write G^ft = GuGiqG 1 = G i0 GG l = ftG. □ 



Lemma 5.6 Lei Wj(z) and Mj(z) be as in Section \5.1\ It holds that 

wf ) {z)Q j {z)=w^&{z) (5.7) 

^(^(*) = ^yt«,M( a ) (5.8) 

Mf'COg,-^) = Af (fe) (z)G(z) J . (5.9) 
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Moreover, for k > it holds that 



M (9-i+fc) (z) 



1 



-M 



(z)G k+j+1 . 



(5.10) 



(g-l + fc)!^ (g_l)r 

Proof We need the following observation. For a row vector x of appropriate 
length and a scalar y, one has 



because 



(x, y)Gjj — x Gjj-i, 



r — ( 
33 ~ \ 



(5.11) 



We now show (|5.7I) . It is obviously true for j — 0. Assume it holds for some 
j > 0. We get, using (fBTB]) and (fSTTTj) 



= wf ] g 3 G 



Equation (|5 .8[) follows by combining (|5 . T[) and (|4.7[) . whereas (|5.9|) is an imme- 
diate consequence of (|5.7p . Next we compute, using t|5 ,9|) and (|4.11|) . 

r (q-l+fc). 



M; 



'(2) M(«- 1+fc )(z) 



(q-l + k)\ * n ' (q-l + k)\ 



_ M(«- X )(z) 
(9-1)! 

which yields ([5~TU|) . 

Let for j <r — 2 the matrix Kj(z) be given by 



□ 



/ F F' 



F F' 



F' 

j 

F; 



1 r>ti) 



-G 

I -G 



-G 
I ) 



(5.12) 
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Here Kj has q diagonal entries Fj and v — q diagonal entries I = I r . Hence K 3 
has dimensions (yr — qj) x (vr — q(j + 1)). A compact expression of K 3 is as 
follows. Let l q be the last standard basis vector of K 9 , t v - q the first basis vector 
of W J ~ q , and S q the shift matrix of size q x q. Then, similar to (I4.13[) . 



Iq ® Fj +Sq® Fj 



— / f 



(8) J r - S^-, ® G , 



Note that the matrices are empty for j > r — 1, ■^trj^r ( -vl 1 q(z) = (0, . . . , 0) 
and that Gyg is empty for j > r. Hence we define 



K r 







\ 




I 


-G 






I -G 






'■■ —G 




V 


I J 


-l))x(u- 


- q)r, whereas for j 



(5.13) 



( I -G 

I -G 



V 



-G 
I J 



(5.14) 



a matrix of size (v — q)r x (y — g)r. 

In what follows, we need the matrices IO(z) — Ko(z) ■ ■ ■ Ki(z), where the ma- 
trices Ki(z) have been introduced in (|5.12|) . (|5.13|) . (|5.14j) . Then A?(z) is of size 
x (yr — q(i + 1)) for i < r — 2 and of size x (y — q)r for i > r — 1. Note 
that X? is always of full column rank. The next lemma extends Equation (|5.4[) . 
obtained for the case v <q. 

Lemma 5.7 Lei < i < r — 1. For < i < k one has 



= (K\ k (z),K 2 lk (z)) 

where K\ k {z) 6 R^- 4 - 1 )-? and ft? fc (z) £ ]ft«x(r-i-l)(i/-g) are exp i lclt i y given 
by 



(5.15) 



(A- 



K* k (z) 



0! ' ' (g-l)! ; ' 

M( fc +«- l - 1 )(z)G(z) 4+1 M( fe +^ l - 2 )(z)G(z) J+1 s 



D! 
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Fori>k it holds that M k (z)K.' l (z) = 0. 

Proof The case k = has been verified in the proof of Proposition 14.61 Let 
therefore k > 1. To prove that the assertion holds true for i < k, we assume the 
right hand side of formula (15.15[) to be valid for Af k (z)jC' l ^ 1 (z) and proceed by 
induction. To that end we multiply it by Ki and verify the answer. As before, 
we denote the j-th block column of K{ by Kf, for j = 0, . . . ,v — 1. We will 
discern the four cases j = 0, j = 1 . . . , q — 1, j = q and j = q + 1, . . . , v — 1. 

Let j = 0. Then the product Af k K i ^ 1 K° becomes j^M^'^F,. The 
analogue of (|5.2j) . with Mi and Ft substituted for M and F, yields that this 
equals ^ k _]_ 1 y M^ l ^ 1 \ as should be the case. 

Let 1 < j < q — 1. One obtains 

^"^wV^Vt^ + ^S^ (5 - 16) 

The analogue of Equation JO]) yields M/ fc+J ' _<_1) F/ = - M ^l-' . Hence the 
right hand side of (|5.16[) becomes 

11, -7 .^+,-^ +Mf +,-^ 1 lM^ Fi 



(k-iy.j\ K k+j-i 1 1 ' (k-i-l)lj\ k + j 



(fc + j-i-l) 



Invoking the analog of (|5.2I) again, we can rewrite this as ^ k _j_ 1 y ~ iJ ~ L j} j a 

typical block of ,H\ k , as required. 

Next we consider the more involved case j = q. In this case the block column 
K\ has entry —QiG on the (q — l)st row (see Lemma l53j) and I on the g-th row. 
Hence we get 

1 JI / f( fc -*+9- 1 ) juf{k-i+q) 

Ar ^' i5 ■ = (^7ji<-i^TT)r ^; • G+ ^^ G •)■ < 5 - 17 > 

Using dSHJ) we obtain M^ l+q ^ 1) g i G = M^ i+q ~^> G i+1 . In view of gjTJ, it 
holds that M^ i+q) = (q + k- i)M^ i+q ~^G. Hence we van rewrite the right 
hand side of (|5.17[) as 

1 , M (k-i+q-l) M (A,_i+,_l) 

( - —7 TTi + (9 + k - l) : )G l+1 



(fc-i)! v (9-1)! ^ ' 9 ! 

which is equal to 

1 M^+g- 1 ) ^ 



the first block of 72.?, , as was to be shown 



Finally we treat the case q+l < j < v — 1. The block columns if| have — G 
at the (j — l)st row and I at the jth row. Hence, we obtain 

N 1C K{ = -i— - M .. . G + ., G ). (5.18) 

(fc-i)! v (j-1)! j! ' 
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Since k-i+j-1 > q, we apply fiTT]) to get M^~ i+ ^ = (k-i+j)M^ k - i+ ^-^G, 
and the right hand side of f|5 . 18[) reduces to 

1 M^ k ' i+j ^G i+1 



(k-i-iy. j\ 

a typical block of 1Z? k , as desired. This settles the proof of the validity of 
Equation (f5T5]) . □ 

Theorem 5.8 It holds that Af k (z)iO(z) = 0, for i > k. For fi < r - 1, the 

matrix iC ,1 ^ 1 (z) is of size vr x (vr — qfi) and has full rank, equal to vr — q[i. 
If H > r - A^ Al_1 (z) is of size vr x (y — q)r and has full rank, equal to (v — q)r. 
Summarizing, the kernel of N(z) is (vr — qmm{fj,,r})-dimensional and spanned 
by the columns of JC ,1 ~ 1 (z). The rank of Af(z) is equal to qmin{/j, r} < vr and 
therefore M(z) never has full column rank. 

Proof We show that N k '(z)iC (z) — for i > k, for which it is clearly sufficient 
to show that M k {z)K k {z) — 0. Starting point is Equation (|5.15l) for i = k — 1. 
We have 

Af k (z)iC k - x (z) = 
M { k °\z) M^ 1] (z) M^(z)G(z) k M^^{z)G{z) k , 



( 



0! ' ' (g-1)! ' 5! ' ' (v-l)\ 



We multiply this equation with the block columns Kt and, as above, we discern 
the case j = 0, 1 < j < q — 1, j = q and j = q + 1, . . . , v — 1. 

For j = we get M k (z)K, k - 1 Kl = M ( k 0) F k = uw k F k = 0, whereas for 
1 < j < g — 1 one computes 

mU-i) m U) 
M * t k-i Ri = ^-Fi + ^-F k = 0, 



in view of an analogue of (|5.2p . 
For j — q, we obtain 

(q-l)l q\ 
We can now use Equation (|5 . 10|) and (|4.10j) to get 

M^~ 1] Q k G _ M^-^G k+1 _ M^G k 
(9-1)! ~ _ ~~ qY~ • 

Hence, the right hand side of (|5 . 19[) is zero. 

Next we consider the case q+1 < j < v + 1. We then get, parallel to (|5.18p . 

Afkri-i M(J - 1)G , MU \ c k 
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which is zero, in view of Equation (|5.10[) . 

To show that Af(z)jCi i - 1 (z) = 0, one has to show that M k {z)^' 1 {z) = 0, 
for all k < (i—l, but this has implicitly been shown above. The other statements 
in the theorem have already been addressed before. The theorem is proved. □ 



5.3 The case v > q + r 

We follow the approach leading to Proposition 14. 71 We observe that the matrix 
Af(z) for v > q + r can be decomposed as 

J\f(z) = (W*(z) O^qxrlv-q-r-l)) , 

where Af*(z) is the u Af(z) matrix" for the case v = q + r—1, since all derivatives 
of M(z) of order higher than q + r — 2 vanish. Let be the 

matrix for the case v = q + r — 1 . Put 

fr-i = tf!T l o 

where I is the identity matrix of order r(y — q — r — 1). If (i < r, then K,* (z) 
is of size ((q + r — l)r x (r(r — 1) + (r — (i)q), and if (i > r, then it has size 
(q+r— l)r X r(r— 1). Then JC^ 1 (z) has size vr x (ur — (iq) for /i < r and has size 
vr x {v — q)r for (i> r. In short, KT' (z) has dimensions 2/r x (i/r — q min{/i, r}). 



Theorem 5.9 Let v > q + r. TTie kernel of the matrix M(z) is spanned by the 
columns o/ A^ -1 (z), has dimension vr — (iq if (i < r and dimension (y — q)r if 
(i>r. So dimkei (J\f(z)) = vr — qmm{(i,r}. 

Proof Similar to the proof of Proposition l4.7l using the results of Theorem 15. 81 
for the case v = q + r — 1. □ 



6 Intermezzo, properties of A°(z) 

The results of this section will be used in Section [71 where we want to find 
(special) right inverses of the matrix A4°(z). 

We focus on the matrix A = (A , . . . ,A U - 1 ) £ R? xi/r , the first block row 
of A, the matrix defined in Section [2j One directly sees that the rank of ^4° is 
equal to min{g, v}, although it also follows from Theorem [23] with (i = 1, Hence 
.A has full rank iff v > q. We introduce the matrix B° £ M. l/rxq consisting of 
the r x q blocks B k as follows. 

/ B° 

B°=l : J Ki.1. ) 
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where each B k has elements 



'•' [ I) Hse. 

for i = 0, ...,r- 1, j = 0,...,g- 1. 

Lemma 6.1 Let v > q. Then A has full row rank and A°B° = I. In other 
words, B° is a right inverse of A . 

Proof We have to compute the ij-elements of T := ^2^=0 A k B k . Using the 
definitions of the matrices A k and B k that only have nonzero entries on corre- 
sponding anti-diagonals, we see that A k B k is a diagonal matrix. Hence we only 
have to consider the iz-entries of T. Note that B k = for k > q. One obtains 

9-1 

T lt =J2(A k B k h 

k=0 
9-1 



5 

k + 1 



g! «4/ g -i-A (-i)fc-i 
— 1 — i)! ^ \ / k + 1 

To compute the latter summation, we write it as 

r l 9-1 / -i \ r l / -, -\ 

= / (l-a^-Vdi 
Jo 

= B(g-i,* + l), 

by definition of the /3- function B(-,-). Using the well-known fact that this can 
be computed in terms of T-functions (B(a, (3) = T(a)T((3)/T(a + /?)) we obtain 

B(q-i,i + l) = ^ -j-^-. 

It follows that T,, = 1. □ 



We need some additional properties. 

Lemma 6.2 It holds that X k := B k Sj + S r B k = 0, if k > q or k < r - 1. For 
i/ie case r < fc < q — 1 (which requires q > r) only the last row of this matrix 
is nonzero. In fact, this row is equal to (— ( k 1 1 )eJ_ r , with the convention 
that ei denotes the i-th basis vector ofW (i = 0, . . . , q — 1). 
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Proof Wc compute the zj-element of X k — B k Sj + S r B k . For i = 0, . . . , r — 1 
and j = 0, . . . , q — 1 it is equal to 

q—1 i — 1 

?— 1 r— 1 

= ZJ B *J +1 l {;=J + i} + 1 {''+l=/}- B »*+l,j 

= ^i,j+l 1 {0<j+l<?-l} + l{0<t+l<r-l}-Si+l,2 
= - B ij + l 1 {0<j< (? -2} + l{0<i<r-2}-B i+lj - 

= (fc + l) 1 {'+J + l = fe} ((- 1 )' 1 {0<j<g-2} + (-l) l+1 l{0<i<r-2}) 
= (u+ 1 ^) 1 {»+J + l = fe}(~ 1 ) I ( 1 {0<J<(J-2} - !{0< l <r-2})- 

Clearly, for i = 0, . . . , r — 2 and j — 0, . . . , q — 2, the last expression in the 
display equals zero, as is the case for i = i — 1 and j — q — 1. We next consider 
the two remaining cases, the first being i < r — 2 and j = q — I. Since we only 
have to consider i = k — j — 1, we get i = k — q, which has to be nonncgative, 
so k > q. But then the binomial coefficient is equal to zero. The remaining 
case is i — r — 1. Then we only have to consider j = k — r, the other values 
of j again give zero. Note that this implies that k > r is necessary to get a 
nonzero outcome, whereas we already know that also k < q — 1 is necessary. 
Hence nonzero elements in the last row of X k can only occur if r < q — 1. 
Under this last condition we find X k _ l rj = (— l) r_1 l{j=fe-r'}- Hence the 

bottom row of X fc equals (— l) r_1 (l{fc =r }, . . . , l{fc=r+ 9 -i}), which is equal 
to ( « )(-iy- 1 el_ r ,{ovk = r ) ...,q-l. ' □ 



Remark 6.3 Here is an example where X as defined in Lemma 16.21 is not 
equal to zero. Take k = r = 1 and q = 2. Then B 1 = (0 1) and X 1 = (1 0). 

Proposition 6.4 Define H k : E -)• M rx « 6y i? fe (» = W r (z)D- 1 B k D- 1 U q (z). 
Then Hk is a constant mapping, Hk(z) = D^ 1 B k D^ 1 , under the condition 
k > q or k < r — 1. 

Proof First we prove the following auxiliary results. One has 

U' q (z) = U^D.SjD; 1 (6.2) 
W' r (z) = D~ 1 S r D r W r (z). (6.3) 

Equation (|6.2p follows from the definition of U q (z) and the elementary identity 
U' q (z) = U q (z)Sj . Equation (|6.3|) can be proved similarly. 
We now compute 

H' k {z) = W;{z)D; 1 B k D- 1 U q (z) + W r {z)D; 1 B k D- 1 U q (z) 

= D- 1 S r D r W r (z)D- 1 B k D q - 1 U q (z) + W r (z)D^ 1 B k D q 1 U q (z)D q S q T D^ 1 , 
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according to Equations (I6.2p and (|6.3[) . Putting S r — D r 1 S r D r and S q — 
DqSqDq 1 , we see that H k satisfies the linear differential equation 

H' k {z) = S r H k (z) + H k (z)Sj. (6.4) 

This equation has a unique solution and we claim that it is given by the constant 
function as asserted. To that end we check 

S r D- l B k D~ l + D~ 1 B k D~ 1 Sj = D~ l {S r B k + B k Sj)D q - 1 

= 0, 

by Lemma 16.21 since k > q or k < r — 1. Furthermore, we have H k (0) = 
D~ 1 B k D~ 1 , since 17,(0) = I q . □ 

Remark 6.5 Equation (|6.4[) has as the general solution 

H k (z) = exp(S r z)iJ fe (0) exp(Sjz), (6.5) 

where the exponentials can be computed as finite sums, since Sj and S r are 
nilpotent. Elementary computations yield for instance that the ij-element of 
exp(Sjz) is given by Q^z 2- - 7 for i > j and zero otherwise. Hence we obtain 

exp(S q z) = U q (z), which is in agreement with the definition of H k (z). 

An example of a solution that is not constant is obtained for r = 1 and 
q = 2. For the case k — 1 one finds directly from the definition of H k (z) that 
•ffi(O) = B 1 = (0 1) and H k (z) = (z l) in view of Remark O This is in 
agreement with Equation (16. 5[) , whose right hand side is equal to 



7 The equation Ai°(z)C = I 

We return to one of our original aims, finding a right inverse of the q x vr 
matrix A^°(z) = (M(z), M^ v ~ l \z)) G R qxri/ . Recall from Theorem O 
and Proposition 13.11 that A4°(z) has rank equal to min{z/, q}. Hence the matrix 
is of full rank iff v > q. Equations like A4°(z)X = b will in general not have a 
solution X for a given b £ IR yrxl , if v < q. In fact, we are interested in solutions 
X that are independent of z. It is easy to see that such solutions only exist if 
6 = and then X — 0. The uninteresting case v < q will therefore be ignored 
and the standing assumption in the remainder of this section is v > q. Under 
this assumption, there are two subcases to discern, r > q and r < q. 

Proposition 7.1 Assume that r > q and v > q. Let I q be the q- dimensional 
unit matrix. There exists a constant (not depending on z) matrix C £ ]g^ rx i 
such that M°(z)C = I q for all z. The equation M°(z)X = b for b £ R q then 
has the constant solution X = Cb. The constant matrix C is unique iff is = q. 
In all cases one can take C = (I v ® D~ 1 )B°D q ~ 1 , with B° as in \6.1\) . 
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Proof In this proof we simply write I for I q . Suppose that we have found a 
constant matrix C with the desired property 



M°(z)C = I. 

By differentiation of (|T. 1 [) k times, with k = 0, 
definition of M.{z) with fi = r, that 



M(z)C = 



/A 





W 



(7.1) 

. , r — 1, we obtain, recall the 
(7.2) 



We note that now M(z) is of size rq x rjv and has rank equal to rq. Hence M. (z) 
has a right inverse, AA(z) + say, and a true inverse in the case that v = q, see 
e.g. Corollarv l3.6l It follows that C should be the first block-column of M.(z) + . 
Next we use the factorization (|3.1[) and note that also A has a right inverse, A + 
say. Then 

X(z)+ = (J„ ® ^(z)- 1 )^/,. ® C/^z)" 1 ). 



Hence, we can choose 



/A 





which means that C is the first block-column of M.(z) + , so 



a = (/ v ®w r (*) _1 )>t + 



o 



V o J 



{I v ®W r {z)- l ){A°) + U q {z)-\ (7.3) 



where (-4°) + is a right inverse of the matrix A°, since A is the first block-row 
of A. But, a right inverse of A is in Proposition [6J] shown to be B°. Therefore, 
we can now explicitly pose our candidate for C, 



C = {I®W r {z)- 1 )B a U q {z)- 1 , 



(7.4) 



where £>° as defined in (|6.ip . Hence we have to show that 

(1) the matrix C in (|7.4|) . in fact doesn't depend on z, 

(2) A / l°(z)C = /. Using matrices introduced in Section [31 we write 

C = (7<g> (Wriz) -1 D~ 1 ))B° D^Uqiz)^ 1 
= (I <g> W r (—z)D~ 1 )B°D q ~ 1 U q (—z). 
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Decomposing C as 



where each block C (fc = 0, . 



C 



/ C° 

1) is of size r x q, we get 
-i d& n— if 



= W r (-z) ' D-'B*D-'U q {~z) 



Hence we see that C = Hk(—z), which was in Proposition 



_ shown to be 

constant and equal to D^ 1 B k D^ 1 , if we have k < r ~ 1 or fc > q. Obviously, 



this is true of k = 0, . . . , r — 1, but for fc = r, . . . , v — 1, we have fc > r > q by 
assumption. This proves claim (1). Since C is constant in z, we can take z = 
in ClU) . 

For the second one we have 

M°(z)C = U q (z)A°(I ® VK r (z))C 

= ® W r («))(/ ® W r r (2r)- 1 )Bt/ fl («)- 1 

- U q {z)A°BU q (z)- 1 
= 1, 

in view of Lemma 16.11 Finally, if v — q, then A4(z) is invertible, which implies 
that C is the unique constant matrix solving Ai°(z)C = I, since in this case 
Equation (|7.2|) has a unique solution. □ 



Remark 7.2 The special choice (A°) + = £>° in the proof of Proposition 17.11 
is rather crucial in finding a right inverse of A4°(z) that doesn't depend on 
z. We illustrate this with the following example. Our point of departure is 
Equation (|7.3[) with jj, = l. 

Recalling (|2~T0)) . we can take (A°)+ = £ v , r (0) _1 .A + £i, 9 (0) - \ with .4+ any 
right inverse of A We choose „4 + = „4 T and compute 



B := C r (0)-M 1 £i,,(0)- 1 = £„, r (0)-M 1 e 



because £i j9 (0) = Since .Aj = /oej (Theorem I2.3[) and £i/, r (0) j3 - = 

(5)(-Sr) <_i , for the fc-th block B fc of g we get B fc = ([) (-l)'-*/*,,^ = 

f ej (fc = 0, . . . , v - 1). We conclude that 6 = .4 T . 

In order to see that this may result in a right inverse of A4°(z) that depends 
on z, we choose q = r = v = 2. The conditions of Proposition 17.11 are then 
satisfied. We have 

z 1 
z 2 1 2z, 



and it follows from the above that 



/ 1 


°\ 








—z 


1 
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We close this remark by noting that there also other right inverses of M°(z), 
depending on z, but still having a simple structure. An example (essentially 
taken from [5], where it was only given for v = q = r in a slightly different 
situation) is 



\ fJ (l'-q)xq/ 



where fo is the first basis vector of W. This follows from the easy to verify 
identity 

M°(z)(I v ®f ) = (U 9 (z) (W_ 9) ). 

The assertion of Proposition lTTTl is not true if r < q (the second subcase). Indeed, 
in the proof of this proposition we used the fact that all C k are indeed constant 
matrices, under the condition r > q. If this is not the case, r < q, the matrices 
C k for k = r, . . . , q — 1 are not constant, in view of Proposition 16.41 Let us give 
an example to illustrate this. Consider the case q = v = 2 and r — 1. Then 



and the equation M (z)C = 1% has the unique on z depending solution C = 
C(z)=M°(z)- 1 . 

We now treat the case v > q in more detail. To that end we need the following 
auxiliary result. 

Lemma 7.3 The subspace of the kernel of A4°(z) that consists of vectors that 
are constant in z, i.e. the intersection f] z ker(A^°(z)), is (y — q) + r -dimensional. 
This subspace is equal to the kernel of M.(z) with [i — r, which is the same for 
all z and hence can be parametrized free of z. 

Proof The first observation is that a vector x in ker A4°(z) that doesn't depend 
on z also satisfies A4(z)x = for arbitrary [i^ in particular for fj, = r. The case 
v < q follows from Theorem l5.2[ since in this case the kernel of Ai(z) is the null 
space for all \i > r. 

Let then v > q. Let x be a column vector consisting of r-dimensional sub- 
vectors xo, . . . ,x„_i that don't depend on z. Recalling that A4° consists of a 
row of blocks (uwY n \ we have 

v-l 

M°X = J2( UW ) (n)x n 
n=0 



v — \ n 



u ik) w (n-k) 



EE 

n=0 fc=0 



^ ^ \k 

k—0 n—k 
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The vectors are zero for k > q and otherwise linear independent. Hence, to 
have the above sum equal to zero is equivalent to 



v-l 

E 

n—k 



for k = 0,...,q- 1. 



The equation for arbitrary 1 < k < q — 1 can be differentiated to get 



E 

n—k 



W 



(n+l-fc) 



Xn — Oi 



which, subtracting from the equation for k — 1 yields 

valid for fc = 1, . . . , q. The only constant solutions to these equations are the 
zero vectors, so xq = ■ ■ ■ = x q -2 = 0. 

On the other hand, for k = q — 1 we keep the equation 



E 

m=0 



m + q — 1 
5-1 



Now we relabel the unknowns by setting y m = 1 )a; m + 9 -i to get 

J2 w{m) y™ = °- 



m=0 



We differentiate this equation j times, with j = 0, . . . , r — 1 to get 
/ w (0) ... w (y-q) \ 



\w 



(r-l) 



(v— g+r— 1) 



?y = o, 



(7.5) 



where y is obtained by stacking the y m . The hrst block-column in the above 
matrix is W, the second can be written as S r W, up to the last one equal to 
S"~ q W . Hence the above system of equations can be compactly written as 

(W S r W ■■■ S^- q W)y = 0. 

Let A be the diagonal matrix with elements An = i. A simple computation 
shows that S r W = WS r A, and therefore S^W = W(S r A) k . Writing S r A = T, 
we can rewrite the last equation in y as 



(W WT ••• WT"-i)y = 0. 
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Since W — W r (z) is invertible for any z, this reduces to 



(i T 



T v -i) y = 0. 



Since the coefficient matrix has full row rank equal to r, its kernel has dimension 
(v — q + l)r — r = (y — q)r. Actually, this kernel is spanned by the columns of 
the (v — q + l)r X (y — q)r matrix 



This proves the claim. 



(-T 

I r -T 

Ir 



\ 



-T 

Ir J 



□ 



Remark 7.4 The result of Lemma l7.3l is also valid for /i > r. This can be seen 
from Equation (|7.5[) . Indeed, if /i > r one has to extend the cocfhcicnt matrix 
with additional rows, all involving derivatives w^k), with k > r. But these are 
all equal to zero. 

One might think that the assertion of the lemma can alternatively be proven 
by explicitly computing the matrix K, fJ '~ 1 (z) for \i = r (noting that (D" 1 <g) 
represents the kernel of A4(z) in view of (|1.2|l ) and showing that 
it is not depending on z. It turns out that this idea is false, as shown by the 
following simple example. 

Let q = l, r = fi = v = 3. We compute K.'- l ~ 1 (z) and show that it is not 
free of z. According to the results of Section [5^1 we find for Ko, K\ and K2 the 
following. 

\ 



K (z) 



(- z 








-1 


— z 








Q 


1 


—z 








-1 














1 


























f 











-1 


— z 











1 











-1 














1 


























f 


























f 





V 




















1 



Ki{z) = 



(-* 
1 

if 



0" 

V 
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and the product K?{z) 
by the columns of 



Kv{z)K\{z)Ki(z) yields the kernel of M.(z) spanned 



(D^<S>h)iC 2 (z) 




We see that the last column of fc 2 (z) has a term — 3z in the fourth row, hence 
this parametrization of ker(.M(z)) is not the one we are looking for. The 
reparametrization of kei(M(z)) given by {D^ 1 ® I 3 )K, 2 (z)R(z) with 



R(z) 



(I 
1 
1 

o"o"o 



Vo o o 



5 6z\ 
4 

2"'6"0 
2 
2 / 



yields (Z^ 1 ® h)K 2 {z)R{z) =: t, with 



K = 




which is free of z. The procedure in the proof of Lemma 17.31 (in this case 
y m = x m , since q — 1) yields a kernel of Ai{z) spanned by column vectors that 
are not depending on z. The result of that procedure is the matrix K, above. 

Proposition 7.5 Let q < v < r . Then the equation A4°(z)C = I q admits a 
constant solution C £ W rxq . The dimension of the affine space of constant 
solutions is equal to (y — q)rq. 



Proof According to Proposition 17.11 a constant solution C exists. Any other 
constant right inverse C is such that the q columns of C — C belongs to the 
kernel of A4° and hence to the kernel of A4 for [i — r. In view of Lemma [731 a 
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basis of this kernel can be obtained can be obtained by choosing (v — q)r linearly 
independent vectors. Applying this result to each the columns of C — C, we 
obtain the result. □ 



To illustrate the fact that the constant solution C of Proposition l7.5l is in general 
not unique, we consider the case q = 1, r — v — 2. Then A4°(z) = (l z l) 

and all constant solutions are given by C — C a ,b — (a b 1 — a) T with 
a,i)£l, which form an affine space of dimension (y — q)qr = 2. 

Of course all right inverses of Af°(z), also those that depend on z, are given by 
a much larger affine subspace. Assume v > q and let Cq(z) be any right inverse 
of Af°(z). Then any matrix C[z) = C {z) + X{z), with X[z) G M urxq a matrix 
whose columns belong to keiAf°(z) is a right inverse. Since, dimkerA/"°(z) = 
vr — q, the affine subspace of these right inverses has dimension (vr — q)q. 

The natural extension of the equation M°(z)C = I q is M(z)C = I m , with 
M(z) of order fiq x vr and l m the identity matrix of order ixq. The matrix 
M.(z) has rank equal to min{/x, r} x min-jV, q} and therefore has full row rank 
if and only if v > q and (J, < r. Hence, under the latter condition, and only 
then, a right inverse exists, and the equation A4(z)C = l m has a solution. This 
equation can be decomposed as 

°\ 



ij 

where every Cj has size vr x q. Parallelling our previous aim, also here one 
could be interested in finding solutions C that are constant in z. For Co we 
are in the previous situation, since a constant Co satisfying M°(z)Co = I q , 
also satisfies M. k (z)Co = for all k > 1. The situation for the the other Ck 
is different. Consider for example C\. It should satisfy Ai°(z)Ci = and 
A4 1 (z)Ci = Iq. However, this is impossible for a C\ that is constant in z, 
since differentiating Ai°(z)Ci — yields Ai 1 (z)Ci = 0. We conclude that the 
equation M(z)C — I m for [i > 2 has no constant solutions. 

Nonconstant solutions are for instance Moore-Penrose inverses. These can 
be obtained by using the Moore-Penrose inverse of the matrix A. It follows from 
the proof of Theorem 12.31 that for fi < r and v > q, the matrix A(A) T is the 
identity matrix. Hence (»4) T is a right inverse of A. Using then Theorem 12.21 
and Proposition 13.11 we obtain that 

M(z)+ = {i v <g> Wrizy^c-IA^ c-} q (i, t ® Uqiz)- 1 

is a right inverse of M{z). The inverses Uq(z)^ 1 and W r (z)~ 1 can be computed 
easily, since one has for instance U q (z) — U q {z)D^ and U q (z)^ 1 = U q {—z). 
The inverses C~l and £~* q can be computed in view of the formulas just above 
Theorem 12.21 Since C~l = /^^^(O) -1 , one obtains that its ij-block (i > j) is 
given by Q)(S'J) l_J (— l) 8- -'- Summarizing, we have 



M°(z) 



K M ll - 1 {z) J 



(Co 



Cu-l) = 
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Proposition 7.6 The matrix M(z) has a right inverse iff v > q and [i < r, in 

which case a right inverse is 

M(z)+ = (i„ ® w r (^)- 1 )r I ,, r (o)- 1 ^ T J c M ,,(o)- 1 (^ ® u q {zy\ 

All right inverses form an affine space of dimension (yr — fiq)^q. 
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